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Abstract: The requirement of Yukawa coupling unification highly constrains the SUSY 
parameter space. In several SUSY breaking scenarios it is hard to reconcile Yukawa coupling 
unification with experimental constraints from B(b — > 57) and the muon anomalous magnetic 
moment a«. We show that b — r or even t — b — r Yukawa unification can be satisfied 
simultaneously with b — > 57 and in the non-universal gaugino mediation scenario. Non- 
universal gaugino masses naturally appear in higher dimensional grand unified models in 
which gauge symmetry is broken by orbifold compactification. Relations between SUSY 
contributions to fermion masses, b — > 57 and which are typical for models with universal 
gaugino masses are relaxed. Consequently, these phenomenological constraints can be satisfied 
simultaneously with a relatively light SUSY spectrum, compared to models with universal 
gaugino masses. 
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1. Introduction 



Supersymmetric grand unified theories [SUSY GUTs] are well motivated possibilities for 
physics beyond the standard model [SM] . They provide an explanation of charge assignments 
of quarks and leptons [jl|, ^, J3|. Furthermore the unification of gauge couplings in the super- 
symmetric version of the standard model supports the idea of an underlying theory which 
unifies three seemingly unrelated gauge symmetries at a scale Mq ~2x 10 16 GeV § |,|. 

There is also a possibility that Yukawa couplings of quarks and leptons unify in a similar 
way 0. This is certainly not a necessity since GUT symmetry breaking effects can be 
incorporated into Yukawa matrices at the GUT scale. Furthermore, each of the MSSM 
Higgs doublets may originate from more than one unified Higgs representation in which case 
the relations between Yukawa couplings can be basically arbitrary. (Although given by the 
underlying theory, these relations might be hard to recover from measurements conducted 
at low energies.) Nevertheless the hope is that the underlying theory is quite elegant and 
certain simple relations between Yukawa couplings hold. The simplest and the best motivated 
relation is the third generation Yukawa coupling unification at the GUT scale: \ = X T in 
SU(5) or X t = X b = A r in SO(10). 

A precise analysis of gauge and Yukawa coupling unification requires two loop renormal- 
ization group running and one loop weak scale SUSY threshold corrections. While gauge 
coupling unification is not very sensitive to the exact form of SUSY spectrum, the fate of 
Yukawa coupling unification crucially depends on SUSY threshold corrections H. This is 
due to the fact that gluino and chargino corrections to the bottom quark mass are enhanced 
by tan/3 and can naturally be as large as 50%. Therefore the success of Yukawa coupling 
unification strongly depends on the SUSY breaking scenario under consideration. Alterna- 
tively, looking at it from the other side, requiring Yukawa coupling unification can point to a 
preferred SUSY breaking scenario or to a region of the parameter space within each scenario. 

Similarly b — > s~f and the anomalous magnetic moment of the muon receive contribu- 
tions from SUSY loops. The size of SUSY contributions and their signs depend on the SUSY 
spectrum. Whether SUSY contributions enhance or suppress these observables compared to 
the standard model predictions mostly depends on relative signs of the gaugino masses Mi, 
Mi, M3 and the \i term. 

For example, in the mSUGRA scenario where all gauginos have the same mass M1/2 at 
Mqut, Yukawa coupling unification can be satisfied with negative SUSY threshold corrections 
to the bottom quark mass. The gluino-sbottom loop is the dominant contribution and is 
negative for fxM^ < 0. On the other hand, the chargino contribution to b — > S7 can interfere 
constructively or destructively with the SM and the charged Higgs contribution. Since the 
SM and charged Higgs contribution is already somewhat too large, the chargino contribution 
is expected to lower the branching fraction B(b — > 57). This happens for \xAt < 0, where At is 
the top trilinear coupling. The low energy value of At is related to — M3 due to renormalization 
group [RG] running. Therefore, in order to accommodate the b — > S7 data, fiM^ > is 
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preferred. 1 Finally, the SUSY contribution to the anomalous magnetic moment of muon, 
a M , is preferred to be positive to comply with the data |J. The chargino-sneutrino loop is 
dominant in this case and this contribution is positive for [xM<i > 0. In conclusion, in the 
mSUGRA scenario both b — > 57 and seem to prefer fj,Min > while Yukawa coupling 
unification strongly prefer fiM\/2 < 0. This discouraging result led to consideration of Yukawa 



coupling unification within other well motivated frameworks [10, |Tl|, 13, 14, |15| ]. 

In Refs. Jl0], |ll|, [l2| SO(10) Yukawa unification was considered together with SO(10) mo- 
tivated boundary conditions for soft SUSY breaking parameters. It was found that t — b — r 
Yukawa unification can be satisfied with positive fi in a special region where the condition 
Aq = 4mf 6 = 2mf between universal trilinear coupling (Aq), universal squark and slep- 
ton masses (m\§) and universal Higgs masses (mio) is approximately satisfied. Additional 
splitting of Higgs masses (rriH u , mH d ) at the level of 10% is necessary in order to obtain 
electroweak symmetry breaking radiatively in large tan/3 regime and large m\Q (> 1 TeV). 
The reason for Yukawa coupling unification to work in this region is that chargino corrections 
to the bottom quark mass are enhanced and can dominate the gluino correction, leading to 
small or even negative total SUSY threshold correction to the bottom quark mass [1C]. This 



region also has other very compelling features. The masses of first two generation scalars are 
large, the order of uiiq, which suppresses flavor and CP violation (and also proton decay), 
while the masses of the third generation squarks and sleptons are below 1 TeV keeping this 
region natural with respect to electroweak symmetry breaking. Although this looks like a 
drastic departure from mSUGRA scenario, it actually may originate from it when additional 
RG running above the GUT scale and GUT scale threshold corrections are properly taken 



into account [10] 



In anomaly mediation the gaugino masses M2 and M3 have opposite signs and so it might 
be possible to simultaneously achieve negative SUSY threshold correction to the bottom quark 
mass and positive SUSY contribution to cu. However in this case the chargino contribution 
to b — > sj constructively add to the SM and charged Higgs contribution and thus have to 
be very suppressed. This is possible only if scalar masses are very heavy (at least few TeV), 



which is problematic with respect to naturalness constraints [13|. 

Another approach to accommodate both b— r Yukawa coupling unification and constraints 
from b — > 57 and was taken up in Ref. H] where non-universalities in gaugino masses in 
supergravity models were considered. It was found that Yukawa unification can be achieved 
with an accuracy of a few percent in significant regions of SUSY parameter space. The 
corresponding sparticle spectrum is relatively light, consistent with naturalness constraints. 

Non-universalities in gaugino masses can be very easily obtained (and are quite generic) 
in higher dimensional GUT models in which GUT symmetry breaking is achieved by orbifold 



1 Constructive interference of the chargino contribution with the SM and charged Higgs contribution is not 
ruled out. However in this case this contribution has to be very small which requires very heavy superpartners. 
This is especially true when considering Yukawa unification, since chargino contribution is strongly enhanced 
in large tan f3 regime. Note however, for the same reason even the case with preferred sign of fiAt is strongly 
constrained by b — » sj. 
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compactification of extra dimensions. The doublet triplet splitting of Higgs fields is also 
achieved in an elegant way and proton decay due to dimension 5 operators (which is a serious 
problem in 4-dimensional GUT models |16[| ), can be naturally suppressed in these models. The 
common feature of these models is the existence of a brane or several branes at orbifold fixed 
points on which the gauge symmetry is restricted to be a subgroup of the GUT symmetry. 
This results in an effective 4-dimensional theory with gauge symmetry given by an intersection 



of gauge symmetries on orbifold fixed points |17L 18, IS, 20, [21]. 

The existence of a brane with restricted gauge symmetry plays an important role in 
gaugino mediation. If SUSY is broken on a brane with restricted gauge symmetry, non- 
universal gaugino masses are generated. For example, if using proper boundary conditions 
SU(5) is broken on a brane down to the SM, non-universal gaugino masses Mi, M2, M% can 



be generated on this brane [18|. Even more interesting is the situation for SO(W) models 



in higher dimensions |20|, |2l| which can contain branes with gauge symmetries being differ- 
ent subgroups of 50(10). Using proper boundary conditions, fixed branes with Pati-Salam 
SU(4) x SU(2) L x SU(2) R , Georgi-Glashow SU{5) x 17(1), or flipped SU(5)' x U{1)' gauge 
symmetries can be obtained. If gauginos get masses on these branes, the gauge symmetry 
relates gaugino masses of the MSSM at the compactification scale. In the case of Pati-Salam 
gauge symmetry, gaugino masses M2 and M3 are free parameters while the M\ is given by a 
linear combination of M2 and M3 [^]. The case of Georgi-Glashow gauge symmetry leads to 
universal gaugino masses and finally in the case of flipped SU(5)' x U(l)' gauge symmetry, 



M2 = -M3 and Mi is an independent parameter [22]. If matter fields are localized on a brane 
with GUT symmetry Yukawa coupling unification is expected as in four dimensional models. 

The compactification scale M c in these models is below (but close to) the GUT scale and 
the boundary conditions with negligible sfermion masses and trilinear couplings are realized 
at this scale. Scalar masses and trilinear couplings receive large contributions from gaugino 
masses through the renormalization group (RG) running between M c and the electroweak 
(EW) scale. These contributions are flavor blind and therefore the resulting soft SUSY 
breaking terms at the EW scale cause only a modest flavor violation originating from the 
Yukawa couplings. 2 



In the original works on gaugino mediation [24, 25], the compactification scale was as- 
sumed to be at or above the GUT scale, in order to preserve the success of gauge coupling 
unification. Therefore all gaugino masses are equal at M c . For M c = Mqut, however, this 
scenario predicts the lightest stau f% to be the lightest supersymmetric particle (LSP), which 
violates cosmological bounds on the existence of stable charged or colored relics from the Big 
Bang in models which conserve i2-parity. The situation is different for M c > Mqut- In this 



case additional RG evolution takes place between M c and Mqut |25|- This running generates 
non-vanishing scalar masses and trilinear couplings at the GUT scale. Most importantly, the 
stau mass receives a positive contribution which eventually can make the r\ heavier than the 



2 For discussion of other SUSY breaking scenarios in higher dimensional GUT models, see for example 



Ref. 



- 3 - 



lightest neutralino Z\. This removes the unpleasant charged LSP feature of the scenario with 



M c = M GUT |25 



This cure of the stau LSP problem doesn't apply to the case of non-universal case (in 
which M c < Mqut)- However, in this case the non-universal gaugino masses help us to obtain 
viable SUSY spectra with a neutralino LSP, at least in some regions of model parameter space. 
A recent study |27| of the non-universal gaugino mediation scenario delineates the allowed 
regions of the SUSY parameter space consistent with neutralino LSP and constraints from 
b — > 57 and for various boundary conditions on gaugino masses. 

In this paper we study to which extent Yukawa coupling unification can be satisfied 
together with phenomenological constraints from b — > 57 and in the non-universal gaugino 
mediation scenario. In Sec. ||] we present basic results of non-universal gaugino mediation 
and obtain approximate formulas for gaugino, squark and slepton masses. In Sec. |3] we 
study SUSY threshold corrections to the bottom quark mass based on the SUSY spectrum. 
Understanding of threshold corrections helps us to understand numerical results presented in 
bee. I The summary of our results and our conclusions are given in Sec. [|. 

2. SUSY spectrum of non-universal gaugino mediation 

Gaugino mediated SUSY breaking is quite economical. In the non-universal scenario this 
mechanism is parametrized by three soft SUSY breaking gaugino masses at the compactifi- 
cation scale M c : M a (M c ), a = 1,2,3. Depending on the localization of the Higgs fields, soft 
SUSY breaking Higgs masses and the /U-term can also be generated. Soft SUSY breaking 
masses of squarks and sleptons and trilinear couplings are negligible at M c . More details on 
possible models are presented in Ref. |27| and references therein. 

It is an useful exercise to obtain approximate analytic formulas of the SUSY spectrum. 
These will help us in understanding of SUSY threshold corrections to the bottom quark mass 
and thus the region of SUSY parameter space where Yukawa coupling unification can be 
satisfied. The complete set of two loop MSSM RG evolution equations [RGEs] can be find 



in [28|. In this section, we use one loop RGEs for the gauge couplings, the gaugino masses 
and the scalar masses. We neglect the contribution of scalar masses and trilinear couplings 
in the running. 

From the one loop RGEs for gaugino masses and gauge couplings, 

tlM " 2C a g 2 a M a , C a = _L(33/5,l,-3), a = 1,2, 3 (2.1) 



= C a gl (2.2) 

a log \i 

we obtain the well known result that gaugino masses scale with the square of gauge couplings: 

0, \ = const = 5 , (2.3) 
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where g is the gauge coupling constant at M c . Using weak scale values: a^(M z ) = 0.118, 
a em (M z ) = 1/128, sm 2 w = 0.22, a 2 = a em / sin 2 9 W ~ 0.035, a x = (5/3)a em /cos 2 9 w ~ 



0.017, and the GUT scale value of the gauge coupling, etc 
weak scale values of gaugino masses: 



g /Att ~ 0.04, it is easy to find 



Mi(M z ) ~ 0.4 Mi, 
M 2 {M Z ) ~ 0.9 M 2 , 
M 3 (M Z ) ~ 3.0 M 3 . 

Similarly, one loop RGEs of squarks and sleptons are given by 



dm 2 - 
dlogfj, 



16vr 2 



(2.4) 
(2.5) 
(2.6) 



(2.7) 
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where dots represent terms which are zero at M c , being proportional to scalar masses and 
trilinear couplings. For comparable values of gaugino masses we can neglect contributions 
from terms proportional to Mi (except in the case of right handed sleptons) since these are 
suppressed compared to terms proportional to M 2 and M 3 by smaller value of gauge coupling 
gi and also the group theoretical factors. 

In this approximation the solution can be written as: 



m 2 D {n) 
m 2 L (fi) 
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1 o 
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= -Mi 
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(2.8) 
(2.9) 
(2.10) 
(2.11) 
(2.12) 
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and for the weak scale values we obtain: 



m,Q{Mz 



2.6 Ma 



mu(M z ) = m D {M z ) 
m L (M z ) = 0.6 M 2 , 
m E (M z ) = 0.39 Mi. 



2.6 M 3 , 



(2.13) 

(2.14) 
(2.15) 
(2.16) 



We also need to estimate the top trilinear coupling which enters the chargino contribution 
to the bottom quark mass. The one loop RGE is given by 



1 



■PA t , 



where 



Pa* 



dA t _ 

dlog/i 16-7T 



32 26 

—g 2 M 3 + 6g 2 2 M 2 + -^9iM x + ... . 



In the same approximation as for the squarks and sleptons, we obtain 



16 

A t (jj) = —M 3 



and the weak scale value is given by 



93 (m) 



3M 2 



92(p) 



A t (M z ) = -3.5 M 3 - 0.4M 2 . 



(2.17) 
(2.18) 

(2.19) 
(2.20) 



These results for gaugino, squark and slepton masses should be interpreted with care. 
However the approximation used is precise enough to point out the main features of the 
SUSY spectrum in non-universal gaugino mediation. First of all, from Eqs. ( |2.4[ ) and ( |2 . 1 6| ) 
we immediately see the well-known problem associated with gaugino mediation, the stau being 
LSP. Neither the bino nor the right-handed stau are mass eigenstates however. The lightest 
neutralino is a mixture of gauginos and Higgsinos, and for f\ the D-terms as well as left-right 
mixing have to be properly taken into account. In result, there is quite significant region of 
parameter space with neutralino LSP, as illustrated by Fig. [|. The allowed region opens up 
with increasing M\ , either because of the large portion of Higgsino in the lightest neutralino 
(along the line where electroweak symmetry breaking is possible, where the /i term is small) 
or the large portion of wino in the lightest neutralino (for small values of M 2 compared to 
Mi ) . On the other hand with increasing tan /3 the r Yukawa coupling also increases and the 
terms proportional to Yukawa couplings in RG equations become more important. These 
terms suppress the stau mass and so the allowed region is shrinking. For more results and 
discussion see Ref. ]27|. 

The second important observation is that the gluino, stop and sbottom masses are very 
close to each other and approximately equal to 3M3. This will be crucial for discussion of 
SUSY threshold corrections to the bottom quark mass in the following section. 
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Flipped SU(5): M 2 (M C )=M 3 (M C ) with M Hu JMJ = 0,\i>0 
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Figure 1: Typical parameter space slices of flipped SU(5) gaugino mediation models. Areas hatched 
in red are excluded by the presence of tachyonic particles or by the lack of radiative EW symmetry 
breaking. The magenta areas represent regions with staii LSPs. Blue and yellow areas are excluded 
by LEP. 

3. SUSY threshold corrections to the bottom quark mass 

The dominant corrections to the bottom quark mass come from gluino-sbottom and chargino- 
stop diagrams. Pieces proportional to tan (3 can be approximated by: 



Sm b \ 9 2a 3 nig 
~ ~~q — A 1 tan/? TJT^ 



and 



SUSY 



m b 



A,: A t tan P 



16vr 2 



M 2 



(3.1) 



(3.2) 



Here Msusy is the typical mass of particles in the corresponding loop. In the previous section 
we found that 



Msusy ~ m§ ~ m~ b ~ m~ t ~ 3M 3 , 
A t 3.5M 3 - 0.4M 2 . 



(3.3) 
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Taking aj,{Mz) = 0.118 and X t (Mz) ~ 1, we can obtain very simple formulas for the SUSY 
threshold corrections to m b : 



f— V -^rrMtaii/J, (3.4) 
V m b J 100M 3 v 7 



tfmA** 1 / 1M 2 \ 



m 6 / 5OOM3 V. 9 M 3 J 

Comparing the last two expressions we find 

'6m b \ x± 1 { 5m b \ 9 ( 1 M 2 



1 + -—=■ //tan/3. (3.5) 



m b J 5\m b J \ 1 + 9M 3 J- (3 ' 6) 

Clearly, the gluino corrections will dominate unless M 2 ~ 4OM3. So for positive fi, the 
Yukawa unification can be expected only for negative M3, and the dependance on M 2 and Mi 
is not expected to be significant. The results shown in the figures of the next section are in 
agreement with these findings. Although \xMj, < is strongly preferred as in the mSUGRA 
scenario, non-universal gaugino masses can help to satisfy constraints from b — > 57 and in 
some portion of SUSY parameter space. 

4. Yukawa coupling unification 

In this section, we present results of our numerical analysis. There are several programs 
readily available for numerically solving the necessary, coupled system of RGEs at the two- 
loop level. The most popular ones are IS A JET [29], SoftSUSY pC|] , Spheno [31] and Suspect 



[3^1 . Recently, a satisfactory agreement was demonstrated between the latest versions of 
these codes fl3^| . To calculate the sparticle spectrum and the Yukawa couplings between M c 
and the weak scale, we use ISAJET version 7.64. The calculational procedure implemented 
in ISAJET, with special emphasis on Yukawa couplings and the theoretical uncertainties in 
their determination, is described in a recent publication fll2| . Due to the known theoretical 
uncertainties, we treat our numerical results with a reasonable flexibility, and conclude that 
unification takes place when the GUT scale Yukawas agree within a few percent. For the 
evaluation of the branching fraction B(b — > 57) and the muon anomalous magnetic moment 
a^, we use the calculations presented in Refs. |34| and |55), respectively. 

In our plots, we present theoretically, phenomenologically and experimentally acceptable 
models. We discard models with tachyonic particles, no radiative EW symmetry breaking 
[REWSB] or a stau LSP as unacceptable. We further require that the Z boson have negligible 
decay rates to sparticles and also impose the following LEP2 constraints on the weak scale 
masses of the lightest gauginos (Z\, W\), sleptons (ei,/ii,fi) and Higgs boson |plj] : 

> 37 GeV, > 100 GeV, 

m Sl > 92 GeV, m Al > 85 GeV, m fl > 68 GeV, (4.1) 
m h > 91 GeV. 
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Models that do not satisfy the above criteria are excluded from our study. Models passing all 
the above criteria are further analyzed with respect to the indirect experimental constraints 
from B(b — > sj) and a„. We impose the following limits: 

-10 < a M x 10 10 < 60, and 2 < B(b S7) x 10 4 < 5. (4.2) 



These limits correspond to the experimentally allowed regions approximately at 2a [37]. In 
our plots, we color models that do not satisfy this B(b —* S7) (or a„) constraint by shades 
of red (or gray). Models that pass all the above requirements are marked by green on our 
figures. 

As discussed earlier, for gaugino mediation only the gaugino and Higgs masses are non- 
zero at M c . In our numerical analysis we assume M c = 10 16 GeV. The rest of the soft breaking 
parameters, namely all squark and slepton masses and trilinear couplings are assumed to 
vanish at M c . The only other remaining free model parameter is tan /3, since our requirement 
of the REWSB fixes the rest of the Higgs sector. Since we scan through both positive and 
negative values of gaugino masses, the sign of fi is not an independent parameter and we 
present results only for fj, > 0. 3 Thus, the parameter space of the non-universal gaugino 
mediation scenario is spanned by six parameters: 

Mi(M c ), M 2 (M C ), M 3 (M C ), M Hu (M c ), M Hd (M c ) and tan/3. (4.3) 

This corresponds to the case when gaugino masses are generated on a brane with SM gauge 
symmetry. If SUSY breaking happens on a brane with larger gauge symmetry, gaugino masses 
are further restricted [2C]. In the case Pati-Salam symmetry gaugino masses satisfy 



M x = \m 2 + ^M 3 , (4.4) 


and in the case of flipped SU(5) 

M 2 = M 3 . (4.5) 

In what follows we present results for b — r and t — b — r Yukawa coupling unification in the 
SM scenario and the flipped SU(5) scenario. In both cases we also show results for vanishing 
Higgs masses, Mh u (M c ) = MH d (M c ) = 0. Finally, we briefly comment on the Pati-Salam 
scenario. 

4.1 Bottom - tau Yukawa unification with independent gaugino masses 

To quantify the amount of b — r unification, we use the variable 

Sy b -r = 100 - l\ , (4.6) 



3 The MSSM Lagrangian is invariant under the simultaneous sign change of the gaugino masses, the A and 
B parameters and fj,. Thus, the results for fi < are equivalent to those for n > with gaugino masses taken 
with opposite signs, since the trilinear couplings are zero at M c . We thank X. Tata for emphasizing this point. 
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Figure 2: A scan of the full parameter space with independent gaugino and Higgs masses at M c and 



fi > 0. The variable Syb- T , defined in Eq.(h6), is plot versus the model parameters. Black (red) dots 
mark models that deviate from the central value of the (B(b —* sj)) measurement by 2a. Models 
marked by green dots satisfy all our constraints. 



where y^ and y r are the values of the b and r Yukawa couplings at the compactification scale. 
The factor 100 is introduced so that <5yfc_ T measures the amount of percentage deviation from 
perfect unification. 

To find the regions with the best Yukawa unification, namely regions with 5yb- T ~ 0, we 
scan the full parameter space as indicated by Eq.([4.3|), Fig. || shows the parameter ranges and 
the result of such a random scan. Here, we plot 5y^ T versus the free parameters of the model 
for jj, > 0. We observe two branches of models, differentiated by the sign of M3. As expected 
from the inspection of the SUSY threshold corrections in Sec. ||, the branch with negative M3 
unifies the b and r Yukawa couplings while the other does not. Yukawa unification happens 
in relatively wide parameter ranges. But models that simultaneously comply with unification 
and the indirect experimental constraints are confined close to tan/3 ~ 45, Mi(M c ) ~ 2 TeV, 
-2 TeV < M 2 (M C ) < 1 TeV, and M 3 (M C ) < -1 TeV. In contrast, there is no strong preference 
for particular values of Higgs masses, except perhaps for low Mh u values. This latter fact is 
supported by an independent scan with vanishing Higgs masses Mi lu (M c ) = M^ d (M c ) = 0, 
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Figure 3: Same as Fig. @ except with vanishing Higgs masses at M c . 



shown in Fig. |3]. For this case b — r unification occurs in similar ranges as for non-zero Higgs 
masses. 

Figs. || and ||| point to specific regions of the parameter space where Yukawa unification is 
achieved and the constraints on B(b — > 57) and are satisfied simultaneously. To explore this 
parameter region further, we fix tan (3 = 40 and scan the M2 versus M3 planes for fixed values 
of M\. Fig. U shows the results of these scans. White areas are theoretically unacceptable 
and areas colored by shades of red (gray) are excluded by B(b — ► 57) (a^). Green areas are 
allowed by all the constraints. As expected, in the negative M2 and M3 quadrants there are 
large regions with good Yukawa unification. For M% > 2 TeV significant parts of these regions 
are allowed both by b — > 57 and a^. Thus, Fig. || clearly demonstrates that, for the case of 
independent gaugino and vanishing Higgs masses, there exist parts of the parameter space 
where it is possible to reconcile b — r unification with all the considered constraints. 

4.2 Top - bottom - tau Yukawa unification with independent gaugino masses 

To gauge the amount of t — b — r unification, similarly to Eq.(4.6), we define 



%- 6 - r = 100f=^^-l), ( 1.7) 

\mm{y t ,y b ,y T ) 
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Figure 4: Yukawa unification contours in the M3 versus Mi plane for independent gaugino and 
vanishing Higgs masses at M c . We chose tan j3 = 40 and fi > 0. 



where are the values of the t, b and r Yukawa couplings at M c . Notice that unlike Syb- T , 
Syt-b-T never becomes negative. In Fig. ||, we plot 5yt~b- T versus the relevant parameters 
from a scan for \x > in the high tan (3 region where Yukawa unification is expected. Again, 
we find models with good unification satisfying all the constraints. The conclusion is similar 
to that of the b — r unification case. Provided that we allow high enough |Mi(M c )| values, 
unification is achieved at moderate values of M2 and M3. Since the correlation between 
6yt-b-T an d the Higgs masses seems to be weak just like in the case of 5yb- T , we conducted a 
similar scan for Mh u (M c ) = Mh d (M c ) = 0. The results are shown by Fig. |6|. In a considerable 
part of the parameter space, we find well unified models consistent with all other constraints. 

Both Fig. H and Fig. || show that Yukawa unification is achieved for two narrow regions 
of tan/3, namely 45 and 50. These two regions are correlated with the sign of M2. The region 
with tan (3 ~ 45 corresponds to M2 > while the region with tan (3 ~ 50 corresponds to 
M2 < 0. This can be seen on Fig. |7|, where we present contours of t — b — r unification down 
to 5%. This plot illustrates for tan (3 = 45 that with vanishing Higgs masses and M\ ~ 5 — 7 
TeV there are significant regions of the M2-M3 parameter plane where t — b — r unification 
takes place. In particular, Yukawas unify for M2 > 0. A similar plot, which we do not include 
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Figure 5: A scan of the relevant parameter space with independent gaugino and Higgs masses at M c 
and \i > 0. The variable Syt-b-r is plot vs. the model parameters. Black (red) dots mark models that 
deviate from the central value of the {B(b — > sj)) measurement by 2c Models marked by green 
dots satisfy all our constraints. 



here, for tan (3 = 50 reveals t — b — r unification for M2 < 0. 

In order to illustrate typical sparticle spectra for models with good unification, we show 
a few parameters for selected models in Table |l]. The firs two models have non-vanishing 
Higgs masses and the third model has zero Higgs masses at the compactification scale. We 
find that well unified models typically have the lightest gauginos in the few hundred GeV 
mass range, lightest sleptons and Higgses below 1 TeV, and the lightest squarks in the TeV 
range. This sort of spectrum may allow the Tevatron to cover a limited part of the parameter 
space, while the LHC has a chance to produce several types of sparticles. A LC, depending 
on its center of mass energy and the actual model parameters, may produce only the lightest 
Higgs boson or additionally the lightest gauginos and sleptons. 

As Table [l] illustrates, the lightest neutralino is always almost degenerate with the lightest 
chargino, that is it always has a strong wino admixture. As a consequence, in this model the 
neutralino relic abundance is generally much lower than the experimental value of the cold 
dark matter Q, as it was pointed out in [27]. 
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Figure 6: Same as Fig. || except with vanishing Higgs masses at the compactification scale. 



4.3 More constrained cases 

Results for independent gaugino masses indicate that Yukawa unification happens for a wide 
range of M2, which suggests that it can be achieved also in the case of gaugino masses being 
restricted by flipped SU(5) symmetry (cf. Eq. (|4.5[) ). Motivated by this, we re-scan the full 
parameter space with M2 = M3 and find that the results are qualitatively very similar to 
those of independent gaugino masses, both in the b — r and t — b — r cases. For completeness, 
we present results for t — b — r Yukawa unification in Fig. ||. This is a scan with independent 
Higgs masses although these are not showed in the figure since there is no preference for 
particular values of Higgs masse just as in previous cases. As expected, only solutions with 
tan /? ~ 50 are found in this case. We also performed a scan with vanishing Higgs masses 
with similar results. It is interesting to note that even in this case, when the SUSY breaking 
scenario is characterized by just three soft SUSY breaking parameters: M±,M2 = -M3 and 
tan/3, we find a region with good Yukawa unification as well as b — > S7 and a«. In this 
case Yukawa unification is achieved at slightly higher values of gaugino masses compared to 
the scenario with independent gaugino masses: M 3 < 1.5 TeV, Mi ~ 6 TeV. Typical SUSY 
spectra from this region are represented by model 4 in Table [l|. The spectrum in this case 
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Figure 7: Same as Fig. [| with contours of t — b — t Yukawa unification being shown. On this plot 
tan/3 = 45. 



exhibits similar features as spectra for independent gaugino masses. 

Finally, in the case with gaugino masses constrained by Pati-Salam symmetry ( [4.4| ), we 
do not expect Yukawa unification while keeping a relatively light SUSY spectrum. Based on 
results for independent gaugino masses and the nipped SU(5) case we see that \M\\ 3> I-M2I 
is highly preferred. However, this can be only achieved for large values of M3. Therefore 
Yukawa unification in this case might only be satisfied with gluino and squarks not lighter 
than several TeV which makes this scenario phenomenologically less interesting. 



5. Summary and conclusions 

In this work, we analysed Yukawa unification in models with non-universal gaugino mediation 
of SUSY breaking. We assumed that all soft SUSY breaking terms vanish at the compact- 
ification scale M c (which is set slightly below the usual SUSY GUT scale). The exception 
being the gaugino and Higgs masses, which were assumed to be non-zero and independent at 
M c . We also considered the cases with vanishing Higgs masses and degenerate M2 and M3 
values. 
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Flipped SU(5): MJM C )=M JM C ), \l>0 




Figure 8: A scan of the relevant parameter space with flipped SU(5) related gaugino masses at M c 
and /x > 0. The variable Syt-b-r is plot vs. the model parameters. Black (red) dots mark models that 
deviate from the central value of the a M (B(b — > sj)) measurement by 2c Models marked by green 
dots satisfy all our constraints. 



We showed that b — r, and even t — b — r, unification can be satisfied simultaneously with 
experimental constraints on b — > S7 and a«. This typically happens for \M\\ 3> IM2I, \M$\. 
The large values of M\ are needed in order to have a neutralino LSP in a substantial part 
of the parameter space. This is necessary especially in the large tan (5 region where Yukawa 
unification can be achieved. 

Yukawa unification strongly prefers \xMj, < just as the mSUGRA scenario or SO(10) 
motivated models [[H], [HJ |l2|. The advantage of non- universal gaugino mediation is the 
possibility of achieving Yukawa coupling unification with a relatively light spectrum. This, 
on the other hand, provides a sizeable SUSY contribution to the muon anomalous magnetic 
moment (shown by Table [l]) , while the SUSY contribution to in most of other scenarios is 
negligible |]io| , [lj] , [l3| , [l2j l . The relic density of neutralinos from resulting region of SUSY pa- 
rameter space is typically below the expectations from cosmological observations. For recent 
discussion of neutralino relic density in other frameworks for Yukawa coupling unification, see 
Refs. pL||, M ||. 
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model 


1 


2 


3 


4 


Mi(M c ) 


3929.6 


3778.4 


4596.2 


6360.3 


M 2 (M C ) 


481.3 


777.9 


959.5 


-1613.7 


M 3 (M C ) 


-1532.9 


-866.9 


-1989.5 


-1613.7 


M Hu (M c ) 


950.7 


816.1 


0.0 


996.9 


MhAM c ) 


144.9 


58.8 


0.0 


138.9 


tan /3 


45.2 


44.7 


44.4 


51.1 


ft(MauT) 


0.525 


0.516 


0.524 


0.509 


fb(M GUT ) 


0.522 


0.510 


0.507 


0.507 


fr(M GUT ) 


0.525 


0.532 


0.519 


0.488 


$b—T 


-0.6 


-4.1 


-2.2 


4.0 


&t—b—T 


0.6 


4.3 


3.3 


4.3 


m~„ 


3376.6 


1986.8 


4325.5 


3545.0 




2837.5 


1738.1 


3647.0 


3104.6 


my 
mi 


2864.7 


1721.4 


3642.8 


3044.4 


2500.1 


1494.5 


3231.9 


2649.5 


nir 


2484.0 


1454.2 


3216.9 


2585.9 


nizj 


765.6 


845.9 


1036.9 


1532.4 




1460.5 


1402.3 


1688.8 


2343.8 


niu e 


761.4 


842.1 


1033.8 


1530.3 


m fl 


615.9 


680.6 


911.6 


1357.3 


mv T 


622.5 


679.8 


919.7 


1357.3 


VV 1 


434.3 


640.8 


837.9 


1283.6 


niy 


1610.1 


823.0 


2063.1 


1542.8 


niTi 
£\ 


434.1 


640.4 


837.7 


1283.4 


m h 


122.5 


119.6 


123.1 


123.4 


m A 


843.9 


659.3 


628.2 


877.5 


m H + 


849.9 


666.5 


635.9 


883.9 


M 


1734.7 


875.7 


2253.8 


1661.9 


a M x 10 10 


35.99 


48.21 


18.50 


8.04 


B(b sj) x 10 4 


3.78 


3.65 


3.96 


4.33 



Table 1: Representative model parameters, and (s)particle masses (in GeVs) for models with good 
Yukawa unification. 
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